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Let  u a  suppose  that  at  a  telephone  exchange  calls  are 
arriving  at  the  Instants  X^,  Xg,  where 

0  <  <  Xg  <...<  Va  <...<  oo  •  As  usual ,  in  queueing  theory, 

we  assume  that  the  inter-arrival  times  «n  *  Vn  -  /t^va  (n=l ,2 ,rQ»0) 

are  identically  distributed.  Independent,  positive  random  variables 
with  distribution  function 


1) 

P^en  <  x]  =  p(x),  (n 

Let 

oo 

2) 

«*  =  f  xdP(x)  , 

h 

and 

CO  __ 

f(s)  e  f  ©"“dLFfx) 

Jo 

The  Input  is  said  to  be  a  recurrent  prooess.  We  shall  ass  tarn  that 
there  are  infinitely  many  lines  available  and  that,  therefore,  no 
call  Is  ever  lost.  Ordinarily,  It  Is  assumed  that  the  holding  times 
(the  durations  of  the  connections)  are  identically  distributed, 
independent,  random  variables  with  an  exponential  distribution.  The 
American  folklore,  however,  assumes  that  woman  talk  more  than  men. 

In  this  paper,  therefore,  we  shall  consider  the  possibility  of  calls 
requiring  one  of  several  (s)  different  exponential  holding  times. 

Por  simplicity,  let  ue  consider  the  case  where  s  *  2.  Thus, 
when  a  call  arrives ,  it  will  be  either  of  the  first  kind  with  pro¬ 
bability  p,  or  of  the  second  kind  with  probability  q.  The  holding 


2 


times  of  calls  of  the  first  kind  have*  the  distribution  function 

4)  Hj^x)  -  l-e"*x  . 

The  holding  tines  of  calls  of  the  second  kind  have  the  distribution 
function 

5)  H2(x)  *  1-e’^  . 

We  let  (5(t),  i^(t))  be  a  vector  random  variable  denoting 
the  number  of  calls  of  the  first  and  second  kinds  respectively 
present  In  the  system  at  tine  t.  We  define  (5n»^n)  be 

i»e.  ,  (§n,^n)  is  the  state  of  the  system  immediately  before  the 
n-th  arrival.  Thus,  the  system  is  in  state  Ejk  at  time  t  If 

(5(t)  ,|(t))  *  (J,  k)  . 

We  shall  obtain  the  binomial  moments  of  the  limiting  distribu¬ 
tion  of  the  Imbedded  Markov  chain  In  order  to  accomplish 

this,  we  first  give  a  bivariate  extension  of  Jordan's  inversion 
formula. 

Later,  we  consider  the  general  case  where  we  allow  an  arriving 

call  to  have  one  of  the  s  types  of  holding  time  distributions 

8 

with  probability  pA(  ^  p4  ■  l)c  Each  of  the  holding  tins 
distributions  Is  assumed  to  be  exponential,  i.e., 

-Ptx 

(1  ®  1,2,.  .  „  ,s)  o 


6) 


H^x)  «  1-6 


~3- 

In  section  5 ,  we  show  that  the  dual  problem  of  a  Poisson 
Input  with  the  possibility  of  several  different  holding  time 
distributions  is  easily  reduced  to  the  M/Q/oo  esse  discussed 
by  Takacs  £l]  . 

It  Is  well  known  that  the  probability  theory  of  type  II 
particle  counters  is  intimately  connected  with  that  of  the  queue 
with  Infinitely  many  servers.  We  use  our  results  to  derive  the 
mean  time  between  consecutive  registrations  in  a  type  II  counter 
when  the  particles  arrive  according  to  a  recurrent  process  and 
the  durations  of  the  impulses  produced  are  distributed  as  a 
weighted  sum  of  exponential  random  variables. 

$2.  Extensions  of  Jordan’s  Inversion  Formula. 

If  {Pj^  (k  =  0,1,2,  ...}  Is  a  discrete  probability  distribu¬ 
tion,  ther  the  r-th  binomial  moment  la  defined  by 

00  V 

71  Br*S(')Pk  • 

00  k 

If  the  generating  function  U(s)  *  P^z  is  analytic  in  a  circle 

of  radius  l+€  ,  where  €  can  be  an  arbitrarily  small  positive 
number,  then 

si  B  -4-  affrfel 

p  ^  drr  «-l 

and  the  binomial  moments  uniquely  determine  the  distribution  ^Pk\  . 
Jordan's  inversion  formula  expresses  the  Pk(k  ■  0,1,2,  . In  terms 
of  the  Bp(r  *  0,1,2,  ...). 


Theorem  1  (Jordan) 


-4“ 


9) 


If 


(1^  ■  *  o  »!»•••) 


is  a  discrete  bivariate  distribution,  we  define  the  r^Tg-th 
binomial  moment  by 


10) 


If 


B. 


00  oo 

U(z,w)  ■  y~  21  pv  V  2 
E^O  k2*0  *1*2 


V* 


is  analytic  in  the  region  |zj  <  1+6,  |w|  <  l+€,  then 
11) 


r,+rP 

1  J  TI(s.w) 

B_  _  *  y«$  i —  ■2-= - A- -L- 1 

FiTo  r.  xr.%* 


lr2  12 


~  d* 


(1,1) 


The  fundamental  inversion  formula  now  becomes 

r„-k„  /r. 

(-1)  x  *(-i: 

*•2 


12) 


^1^2  r^Vk^  ^J3k, 


(rlN|fr2>lB 

W  1*2/  rlr2 


Before  proceeding  with  the  proof,  I  shall  state  the  general 

multivariate  result  which  is  proved  similarly. 

Let 


Pfclk2--kn  (kl 


0,1,2,...,...  ,1^  «  0,1,2,...) 


be  a  discrete  n-varlate  distribution.  Defining  the  (r^rg. . .rn)-th 
binomial  moment  by 


13) 


B. 


rlr2* 


oo 


*r_ 


HR)  *”R)?: 


kikg,. .  .k^ 


“5 


th«n  the  proper  inversion  formula  is  seen  to  be 

(14) 


<Il-kl,*<r2-k2>*-- ( v^1)b 


Proof  of  the  Bivariate  Case. 


In  order  to  verify  the  assertion  that 


5  Va  ^^-^2  Wvv’Vg  ’ 


we  substitute  the  value  of  B. 
obtaining 


rlr2 


in  the  right  side  of  (15), 


oo  oo  r,  -k,  r0-k~  /r,\  /r, 

16)  3“  T  (-1)  1  ^(-1)  2  2 


J*l=ki  r2«r^ 

00  OO  00  00 


r-j^k^  r2=K2  m>!rl  n"r2 


Using  the  fact  that 


iy 

Oiann,. 


and 


Interchanging  the  sums 

oo  oo 


n-lc„ 
oo  2 


fpc^Q  ^Eg  l^kgS 


6 


yields 

00 


in-*, 

Vr2*k; 


Recall,  howerer ,  that  (-l)n“^=^  *£  n  *  k 


Thus,  the  only  non-zero  tern  in  our  sum  occurs  when  n^kg  and  m=»k^ 
Then  our  sum  reduces  to 


18) 


f3«  The  Ergodlc  Behavior  of  the  Imbedded  Markov  Chain  (^,  ^n)  » 

IT  we  look  at  our  process  at  the  times  just  before  a  call 
arrives,  then  the  randan  variables  81  •  °)  *'r(^cn  " 

form  a  Markov  Chain  because  vie  have  assisted  exponential  holding 
times.  The  transition  probabilities  are  given  by 

19 )  pjk6n  °  PC^n+l,'^i+l^  “  *  (3#k))  » 

Starting  from  the  initial  diatrlbution  ^  j  t  the  dlstribu- 
tiona  can  be  determined  auecesslvelf  from  the  Chapman* 

Kolmogoroff  equations 

p  (n+1)  -  p  p  (n) 

lm  T&  *jWjk 


20) 


■7- 


It  is  known,  however,  that  it  is  much  more  convenient  to  work 
with  the  binomial  moments  of  the  distribution  (In ,^n) ,  i.e. , 

v?  -«a-  *  ■ 


21) 


Lemma.  Let  X  and  Y  be  two  independent  binomial  distributions 
with  parameters  (p^  ,0^)  (pg,^)  respectively.  Let  Pj^^-P  (x*ki ,Y=kgj 

be  their  joint  distribution.  Then 

“i  V.  /*!Vk; 


22)  B. 


« -A  A  ©fi)v,  ■  G)BKV- . 


The  proof  is  purely  computational. 

Theorem  2,  We  have  «  1  (n  **  1,2,  ...)  and 


23)  B 


(n+l )  _ 


Pll2 


«hx«2'*>  +  P\-X,r2* 


Proof.  By  our  lemma. 


*  <!•*>>  V*j 


2U 


because  under  our  conditions  the  call  is  of  the  first  kind,  with 
probability  p,  and  thus 


-Ax 


and 


is  binomial  with  parameters  J+l  and  e 
Jn+1  *8  binomial  with  parameters  k  and  e”^* 


-8 


Likewise,  the  call  la  of  the  second  kind,  with  probability  q> 
and  thus 


*n+l 

and 

fc+l 


Is  binomial  with  parameters  J  and  e“*x 
is  binomial  with  parameters  k+1  and  e”*^  , 


Therefore,  we  have  proved  (24),  Removing  the  conditioning  on 

€>n  yields 


Multiplying  both  sides  of  (25)  by  '  and  sunning  over  all 
relevant  (J,k),  we  obtain 


26)  B. 


(n+1) 


CBrir'nl  ♦ 


(n) 


(n)  ~| 


rlr2 


It  is  clear  from  (26)  that  if  the  limiting  distribution 
Pj^  =  lim  P^j,^)  *  (J  ,k)}  exists  (which  it  does  in  our  case), 
then  the  binomial  moments  of  ^Pj^  satisfy  the  equation 

<t{r^r2v.) 

27  Brlr2  ~  i-^**xA+r2pJ'  ^*^-1  ,r2  +  q^*i*p2"l  * 

with  BQ0  *  1 

Before  presenting  the  solution  to  the  difference  equations  above 
we  introduce  further  notation.  Let  D(r^,r2)  denote  the  set  of  all 
decreasing  paths  from  (r^rg)  to  (0,0),  A  decreasing  path  is  one 
that  always  goes  down  or  to  the  left.  Let 


9 


°<»^)  - 


If  either  m  or  n  is  different  from  xero;  and  let  C(0,0)  *  1. 
For  any  path  ?  in  D(r1#r2),  form 

C(5*)  *  TTcCbi!**^)  > 

where  the  product  is  taken  over  all  points  of  the  path. 


Theorem  3.  The  solution  to  the  system  of  equations  (27)  is 


28) 


B. 


rlr2 


JY72 

P  9 


The  proof  can  be  accomplished  by  double  induction  but  a  glance 


at  the  equation  shows  that  the  B_  given  by  (28)  fit  the 

rlr2 

equation  and  reduce  to  the  results  of  Takacs  if  p  =  1.  In  his  ease 

there  is  only  one  decreasing  path  as  he  dealt  with  a  1-dimensional 

r 

problem  and  B_  «  TTc«  •  If  we  impose  the  initial  conditions 
r  J=s i  J 


B. 


and 


V-TT*. 


we  see  that  (28)  gives  the  unique  solution  to  our  equation. 
As  an  example  we  compute  Bgg  ° 


There  are  6  possible  paths 
from  (2,2)  to  (0,0).  Namely, 

(2.2) -(2 ,l)-(2,0)-(l ,0)- (0,0) 

(2.2)  -(1,2) -(0,2) -(0,1) -(0,0) 

(2.2) -(2,l)-(l,l)-(l,0)-(0,0) 
(2^)-(2A)-(lfl)-(0,l)-(0,0) 

(2.2) -(l  A)  “(1,1) -(1,0) -(0,0) 

(2,2)-(l  ,2)  *(1,1) -(0,1) -(0,0) 


Consequently , 

B22  *  °22C12C02C01  +  °22C12C11C01  +  C22C12C11C10 

+  C22C21C20C3  0  +  C22C21CllCai  +  C22C2lCUC10  * 

The  General  S-Dlmenalonal  Situation. 

We  now  permit  an  incoming  call  to  have  one  of  s  exponential 
holding  time  distributions.  We  deal  with  the  following  model: 

Calls  arrive  at  a  telephone  exchange  at  times  <  r2  <•  • « <  Tn  <  . » » 
The  inter-arrival  times  0^  *  rn  (n  *  1,2,  . .»)  are  mutually 

independent.  Identically  distributed,  random  variables.  Thus, 


29) 

<  *J  «  p<x> 

00 

30) 

•<  =  f  xdP(x)  <  oo 

Jo 

and 

<(s)  ■  f  e^dPjx) 

3D 

Jo 

Each  call  has  probability  p^  of  having  the  i-th  kind  of  holding  time 
distribution 


-11 


H^fx)  *  1-e 


(1  ®*  1 ,  •  •  •  i  s )  • 


Of  course ,  ^  *  1  .  Aa  usual,  we  assume  that  holding 

time  distributions  chosen  are  mutually  Independent  random 
variables  end  Independent  of  the  sequence  of  times  of  arrival 


We  denote  the  state  of  our  system  at  time  t  by  an  3-vector 
[g1(t)  ,  %2(t) ,  0„o  ,  §8(t)3  where  i^t)  Is  the  number  of  calls 
of  the  l«fch  kind  present  In  the  system  at  time  t  .  We  3hall 
determine  the  binomial  moments  of  the  imbedded  Markov  chain  of 
this  process o 

x-t  %n  -  <s*  )  =  fsj;frn  .  0) . i;<Tn  -  0)]. 

The  transition  probabilities  of  the  s -dimensional  chain  are 
given  by 


3)  p<Jx,....j0>0£1."”’V  “  P®n+1=  <ki>"”>V|5n  * 


As  before ,  It  la  easier  to  work  with  binomial  moments 


Proceeding  as  before ,  we  arrive  at  the  following  equation  satisfied 
by  binomial  momenta  of  the  stationary  distribution 


12- 


<n) 


P.  ,  *  lia  P.  v“', 

Dr»>  00  Jl***Jf 


«£l  riV.t) 

^  Bp,p2..,r  *  a  [^1^,-1  ,r2,. . .  ,r  ^2®*!. 


with 


Boo...o  *  1  0 


l»r2"*»r2»**JP6 


Let  D(p1>  ...  #rg)  denote  the  set  of  all  decreasing  paths 
from  (p1>0..,ra)  to  (0,  ...  ,0),  and  let 

C |...  *®g )  “  g 1 '  '  11  o 

Por  any  path  $  In  Dtr^,,,. ,rg)  form 

0(f)  »  TT  Cta^,...  ^aa) , 

where  the  product  Is  taken  over  all  points  (n^>.«.,m8)  of  the 
path  Jp  o  We  now  can  give  an  explicit  formula  for  the  moments 
of  the  stationary  distribution  \  e 

Theorem  lu  The  solution  to  the  system  of  equations  (3$)  f.s 


(36) 


•13- 


^5.  Remark!  on  the  Poijaon  Input  Case. 


Suppose  the  oalla  arrive  according  to  a  Poisson  process. 
Each  call  will  have  probability  p^  of  having  one  of  the  a 
holding  time  distributions 

s 

HjU)  ^  Pi  *  (ML,  ...  ,s). 

The  only  requirement  on  the  holding  time  distributions  is  that 
they  have  finite  mean,  l.e., 

OO  op  . 

-ij  *  |  xdH^(x)  «  j  (l-H1(x)3dx  <  oo  for  all  i  . 


This  problem  is  trivial  because  we  can  consider  the  input  to  be 
composed  of  a  different  Poisson  processes  each  with  parameter 

APl  .  Thus,  we  can  regard  our  system  as  s  different  M/g/oo 
systems.  If  we  let  I*  (fi^(t)  ,  ...  ,§?(t))  denote  the  state 
of  our  system,  then  by  Theorem  11  of  Reference  (lj  (also  Theorem  1 
page  160  of  (2) ) ,  we  have 


37) 

where 


•Tgi  Ji 

p(§(t)  ■  (k^,  ...  ,ka)J  »  e  *  TT  , 

Si  *  Xpj  (l-Hj  (x)}dx 


(1  m  l>t  o  f  ,s) 


<§6.  A  Related  Particle  Counting  Probl«m, 


In  this  Motion  we  analyse  our  problem  from  the  point 

of  view  of  particle  counting.  We  assume  that  particles 

arrive  at  a  Type  II  oounter  at  times  T^,  vn,...  (n*l 

where  the  inter-arrival  times  ©n  »  -  Tn  are  identically 

distributed,  positive,  random  variables  with  dletribution  function 

F(x)„  As  before,  we  let 

oo 

«<  *  V  xdF(x) 

Jo 

and 

<t(a)  -  7*e”8XdF(x) o 
Jo 

At  the  instant  of  its  arrival  a  particle  produces  one  of  s 
kinds  of  Impulses  „  Bach  particle  has  probability  Pj  of  producing 
the  I-th  kind  of  impulse  the  duration  of  which  has  the  distribution 
function 


38) 


H^x)  *  1-e 


(i  *  1 ,.  •  *  ,8  )« 


Of  course ,  «*  1.  Furthermore ,  we  assume  that  the  durations 

of  the  impulse  times  which  are  chosen  are  mutually  independent, 
random  variables  that  are  also  independent  of  the  sequence  of 
times  of  arrival  ^n]  • 


Although  every  arriving  particle  produces  an  impulse  in  the 
counter,  only  those  particles  that  arrive  when  the  oounter  is  free 
are  registered.  The  times  of  arrival  of  the  registered  particles 


-15 


form  a  subsequence  of  the  sequenoe  |Tn^.  The  times 

between  consecutive  registrations,  9^  *  ,  are 

also  identically  distributed.  Independent,  positive  random 

variables  with  distribution  function  R(x).  In  this  section 

we  shall  compute  the  mean  of  R(x). 

In  this  particle  counting  framework,  the  random  variable 

t  *  (§' ,  a*o  ,  )  denotes  the  number  of  impulses  of  each 

n  n*  n 

type  present  in  the  counter  just  before  the  n-th  particle 
arrives.  In  particular,  R0**»0^  *a  that 

the  n-th  arrival  finds  the  counter  free  and  is  registered. 


The  limit 


11m 

n->oo 


p  (n) 
r0»  * • 0 


■<)».. 


0 


exists  and  is  given  by 


39) 


P 


0*«  *0 


s 


where  the  B„  _  are  given  by  (36).  By  applying  Wald’s 
rr .  .rs 

Fundamental  Identity  of  Sequential  Analysis  (see  [2 "}  page  183)* 


we  conclude  that  the  mean  time  between  successive  registrations 

is 


40) 


16 


Before  we  know  what  kind  of  particle  has  arrived ,  the 
duration  of  the  Impulse  produced  by  the  particle  has  the 
distribution 

s  -» i,x 

41)  H(x)  «  1-  >  p,  e  x 

1*1  1 

Therefore,  the  time  between  consecutive  registrations  In  a 
Type  II  counter,  when  the  particles  arrive  according  to  a 
recurrent  process  and  the  impulse  times  are  distributed  as  a 
weighted  sum  of  exponentially  distributed  random  variables , 
and  the  time  between  oonsecutlve  registrations  In  our  process 

b<VYC. 

has  the  same  distribution.  Hence 


JL 


Is  also  the  mean  time  between  successive  registrations  in  this 
second  process.  Sinoe  any  distribution  function  can  be 
approximated  by  one  of  the  form  (41) ,  this  result  may  be  of 
practical  value  . Unfortunately ,  we  have  not  been  able  to  obtain 
the  variance  of  R(x)  by  the  use  of  the  Imbedded  Harkov  chain. 
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